Dynamical corrections to the DFT-LDA electron conductance in nanoscale systems 
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Using time-dependent current-density functional theory, we derive analytically the dynamical 
exchange-correlation correction to the DC conductance of nanoscale junctions. The correction per- 
tains to the conductance calculated in the zero-frequency limit of time-dependent density-functional 
theory within the adiabatic local-density approximation. In particular, we show that in linear re- 
sponse the correction depends non-linearly on the gradient of the electron density; thus, it is more 
pronounced for molecular junctions than for quantum point contacts. We provide specific numerical 
examples to illustrate these findings. 



Recent attempts to apply static density-functional the- 
ory Q (DFT) to electronic transport phenomena in 
nanoscale conductors have met with some success. Typ- 
ical examples are atomic-scale point contacts, where 
the conductance, calculated with DFT within the local- 
density approximation (LDA), is found to be in excel- 
lent agreement with the experimental values When 
applied to molecular junctions, however, similar calcu- 
lations have not been so successful, yielding theoreti- 
cal conductances typically larger than their experimen- 
tal counterparts [3|. These discrepancies have spurred 
research that has led to several suggested explanations. 
One should keep in mind, however, that these quantita- 
tive comparisons often neglect some important aspects 
of the problem. For instance, the experimentally re- 
ported values of single-molecule conductance seem to be 
influenced by the choice of the specific experimental set- 
up [J, |a, la, tD|: different fabrication methods lead to differ- 
ent conductances even for the same type of molecule, with 
recent data appearing to be significantly closer to the 
theoretical predictions than data reported in earlier 
measurements y. In addition, several current-induced 
effects such as forces on ions Q and local heating 
are generally neglected in theoretical calculations. These 
effects may actually generate substantial structural insta- 
bilities leading to atomic geometries different than those 
assumed theoretically. However, irrespective of these is- 
sues, one is naturally led to ask the more general question 
of whether static DFT, within the known approximations 
for the exchange-correlation (xc) functional, neglects fun- 



damental physical information that pertains to a truly 
non-equilibrium problem. In other words, how accurate 
is a static DFT calculation of the conductance within 
the commonly used approximation for the xc functional, 
LDA, compared to a time-dependent many-body calcu- 
lation in the zero-frequency limit? 

In this Letter we provide an answer to this ques- 
tion. Specifically, we seek to analytically determine 
the correction to the conductance calculated within the 
static DFT-LDA approach and illustrate the results with 
specific examples. A few recent attempts were made 
in this direction. For instance, Delaney et al. 
used a configuration-interaction based approach to cal- 
culate currents from many-body wavefunctions. While 
this scheme seems to yield a conductance for a specific 
molecule of the same order of magnitude as found in ear- 
lier experiments on the same system, it relies on strong 
assumptions about the electronic distribution of the par- 
ticle reservoirs Following Gonze et al. Evers et 
al. ^3 suggested that approximating the xc potential of 
the true nonequilibrium problem with its static expres- 
sion is the main source of discrepancy between the ex- 
perimental results and theoretical values. However, these 
authors do not provide analytical expressions to quantify 
their conclusion. 

Our system is the nanojunction illustrated in Fig. ^ 
which contains two bulk electrodes connected by a con- 
striction. In order to understand the dynamical cur- 
rent response, one must formulate the transport problem 
beyond the static approach using time-dependent den- 
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This yields 



FIG. 1: Schematic of a nanoscale junction between electrodes 
with applied bias V. Due to dynamical exchange-correlation 
(xc) efTects there is an xc field which gives rise to an additional 
voltage drop y'*''" compared to the electrostatic potential dif- 
ference. 



sity functional theory 

[ilElE3 (TDDFT). In the low 

frequency limit, the so-called adiabatic local-density ap- 
proximation (ALDA) has often been used to treat time- 
dependent phenomena in inhomogeneous electronic sys- 
tems. However, it is essential to appreciate that the 
dynamical theory includes an additional xc field beyond 
ALDA |l7l| - a field that does not vanish in the DC limit. 
This extra field, when acting on the electrons, induces 
an additional resistance 7?"^^", which is otherwise absent 
in the static DFT calculations or TDDFT calculations 
within the ALDA Our goal is to find an analytic 

expression for this resistance and then estimate its value 
in realistic systems. We will show that the dynamical 
xc field opposes the purely electrostatic field: one needs 
a larger electric field to achieve the same current, im- 
plying that the resistance is, as expected, increased by 
dynamical xc effects. 

We proceed to calculate the xc electric field. This 
quantity was determined by Vignale, Ullrich and Conti 
using time-dependent current-density functional theory 
(TDCDFT) ^2Qj- In their notation, this field is 



ALDA 
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where E^^^^{r, ui) is the ALDA part of the xc contribu- 
tion, (TQ,/3(r, w) is the xc stress tensor, and e is the elec- 
tron charge. Here n = no(r) is the ground-state number 
density of an inhomogeneous electron liquid. We are in- 
terested in the dynamical effects that are related to the 
second term on the RHS of Eq. (^), i.e., 
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(2) 



In the static limit [31, we can transform i?^^" into an 
associated potential by integrating between points a and 
b inside the electrodes lying on the z axis. We take the 
2;-direction to be parallel to the current flow. The end 
points a and b are in regions where the charge density 
does not vary appreciably, i.e., 9z(l/n)|^ ~ (see Fig.^l. 
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Importantly, we include here only the part of the elec- 
tric field that varies in phase with the current (i.e., we 
take the real part of the stress tensor). In general — 
e.g., at finite frequency — the electric field has both an 
in-phase (dissipative) and an out-of-phase (capacitive) 
component, where the latter is related to the shear modu- 
lus of the inhomogeneous electron liquid. Such capacitive 
components play a crucial role in the theory of the dielec- 
tric response of insulators. We ignore them here on the 
basis that they do not contribute to the resistance [2l|. 
The general xc stress tensor in TDCDFT is given 

by 



(4) 



+({n,uj)V -uSa/B, 



where fjin^uj) and C(n,a;) are the frequency-dependent 
viscoelastic coefficients of the electron liquid, while u = 
j/n and j are the velocity field and the particle current 
density, respectively, induced by a small, time-dependent 
potential. 

The viscoelastic coefficients are given by 
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and 
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(5) 
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where ^^d f^cri^) ^^^^ respectively, the longi- 

tudinal and transverse xc kernel of the homogeneous elec- 
tron gas evaluated at the local electron density n = no(r), 
while e"j, is simply 
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In the representative systems that we examine below, 
the derivatives in the transverse directions x and y ac- 
count for only a small fraction of the total dynamical xc 
field and can hence be ignored. We thus obtain 



We then see that 
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where the viscosity rj = hnitj^o Rer/(a;) is a function of n, 
and therefore of z. The real part of Ci^) vanishes in the 
Hmit of 0. Under the same assumptions of negligible 
transverse variation in current density |22| , we can write 



(10) 



enAr 



where / > is the total current (independent of z), and 
Ac is the cross sectional area jl^l- Substituting this into 
the equation for the voltage drop and integrating by parts 
we arrive at 



ydyn _ 
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(11) 



Because rj is positive - a fact that follows immediately 
from the positivity of energy dissipation in the liquid - we 
see that the right hand side of this equation is negative- 
definite: the electrostatic voltage is always opposedhy the 
dynamical xc effect. We identify the quantity on the right 
hand side of Eq. ITTJ with -R'^^'^I, where i?'^^" is the 
additional resistance due to dynamical effects. According 
to TDCDFT, the current that flows in the structure in 
response to an electric potential difference V is given by 



(12) 



where Gs is the conductance calculated in the absence 
of dynamical xc effects (e.g., by means of the techniques 
of Ref. 0). Solving Eq. l(T^ leads immediately to the 
following expression for the total resistance R = V/I: 



R = 



1 

Gl 



where 
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V 7—dz. 



(13) 



(14) 



The dynamical xc term thus increases the resistance. 

This is the central result of our paper. It shows that 
the dynamical effects (beyond ALDA) on the resistance 
depend nonlinearly on the variation of the charge density 
when the latter changes continously from one electrode 
to the other across the junction. In a nanojunction, this 
correction is non-zero only at the junction-electrode in- 
terface where the electron density varies most. Knowing 
the charge density one can then estimate this resistance. 

Let us thus consider two specific examples that have 
attracted much attention, namely the gold point con- 
tact and the molecular junction formed by a benzene- 
dithiolate (BDT) molecule between two bulk electrodes 
(see insets in Fig. ^ and estimate the error made by 
the (A)LDA calculation in determining the resistance. 
In order to make a connection between the microscopic 
features of these junctions and the density in Eq. (|14|l . 
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FIG. 2: Planar (dashed line) and macroscopic (solid line) 
averages of the charge density for a gold point contact (left) 
and a molecular junction (right). 



we model the charge density n = no(r) as a product of 
smoothed step functions in every direction, i.e., 

no(r) = neE{z,d,^e) 

+nc^lix, h, jc)^{y, w, Jc)^{z, d, 7e). (15) 

The smoothed step function is given by 



g(2+//2)/7 



(16) 



where I is the full-width at half-maximum and 7 is the de- 
cay length. Here, il{z, 1, 7) — 9(z, — /, 7) — Q{z, 1, 7) and 
S(z, ^, 7) = 1 — n{z, 1, 7); Q{z, 1, 7) represents the density 
distribution of the junction ric, which smoothly connects 
to the constant bulk density rie of the two electrodes sep- 
arated by a distance d. Finally, h and w represent the 
lateral dimensions of the junction. 

The electron densities are obtained from self-consistent 
static DFT calculations with the xc functional treated 
within the LDA (25|. The (111) gold surface orientation 
is chosen for both the point contact and the molecular 
junction (see schematics in Fig.EJ. 

In Fig 13 we plot the planar and macroscopic averages 
of the self-consistent electron densities for both systems 
as a function of distance from the center of the junction 
along the z-direction. The macroscopic average is then 
fitted to the simple charge model in Eq. (|15|l . The fitted 
density is then substituted in Eq. H14|l to find the correc- 
tion to the resistance. The estimated value of R'^^^ 
for the point contact is ^^0.2 Kfi (the static DFT resis- 
tance is about 12 Kfi (23]), while for the BDT molecule is 
^AOKU (to be compared with the linear response static 
DFT resistance of about 360 Kfi Hil). As expected, J?'^^" 
for BDT is larger than that for the point contact due to 
the larger variation of the average density between the 
bulk electrodes and the molecule. In Fig. |31we plot the 
resistance in Eq. (|14|l as a function of the ratio Ue/ric 
and the decay constant 7, where we fix Ue to the value 
of bulk gold {vs ~ 3). The resistances of the two spe- 
cific examples are indicated by dots in the figure. It is 
clear that the dynamical contributions to the resistance 
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FIG. 3: The resistance due to dynamical effects as calcu- 
lated from Eq. 1141 with the charge density determined from 
DFT-LDA calculations as a function of the main parameters 
discussed in the text. The resistance of a gold point contact 
and a BDT molecular junction are indicated by dots. 

can become substantial when the gradient of the density 
at the electrode-junction interface becomes large. These 
corrections are thus more pronounced in organic/metal 
interfaces than in nanojunctions formed by purely metal- 
lic systems. 

In summary, we have shown that dynamical effects in 
the xc potential contribute an additional resistance on 
top of the static DFT-LDA one. The magnitude of the 
additional resistance, within linear response and the zero- 
frequency limit, depends on the gradient of the charge 
density across the junction. This additional resistance is 
thus larger in molecular junctions than in quantum point 
contacts. 
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